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The standard model satisfies Sakharov's conditions for baryogenesis but the CP violation in the 
KM matrix appears too small to account for the observed asymmetry. In this letter we explore a 
mechanism through which CP violation can be greatly amplified. First CP is spontaneously broken 
through dynamical effects on bubble walls, and the two CP conjugate phases grow through phase 
ordering. Direct competition between macroscopic regions of both phases then amplifies the micro- 
scopic CP violation to a point where one of the two phases predominates. This letter is devoted to 
a demonstration that spontaneous CP violation may indeed occur on propagating bubble walls via 
the formation of a condensate of the longitudinal Z boson. 



A fascinating development in particle theory over the 
last few years has been the realisation that the standard 
model, possibly with a minimal extension of the Higgs 
sector, has the potential to produce the observed matter- 
antimatter asymmetry in the universe. At the heart of 
the mechanism is the remarkable anomaly structure of 
the Weinberg-Salam theory, which ties changes in the 
topology of the gauge and Higgs fields to changes in the 
baryon number of the universe (l) If these baryon num- 
ber violating processes were somehow biased at the elec- 
troweak phase transition, and subsequently turned off, 
an asymmetry would be generated rather naturally. 

Interest in this idea grew when explicit mechanisms 
producing the required bias were proposed, in minimal 
extensions of the standard model involving extra Higgs 
fields (for reviews see Two simple arguments seem 

to exclude the minimal standard model. The CP viola- 
tion due to the Kobayashi-Maskawa matrix is very small, 
suppressed by many powers of light quark masses and 
mixing angles. A naive estimate assuming analyticity in 
fermion masses indicates that CP violation should enter 
in the combination 



dcp = JUi^j^m^ - TO^)A 
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(1) 



where the product is over like charged quarks, and 
J ^ 10^^ is a product of KM angles. Taking the 
scale A to be the temperature T 100 GeV, one finds 
dcp ~ 10~^^ leading to a baryon symmetry well below 
the required value ^ 10^^". Second, the requirement that 
the anomalous processes be suppressed after the phase 
transition imposes an upper bound on the mass of the 
Higgs boson - estimated from the one loop potential ^] 
to be Mh < 35 GeV, in conflict with the LEP bound 
Mh > 60 GeV. But recent lattice studies |^ indicate 
that the electroweak transition is far more strongly first 
order, relaxing the upper bound on Mh, perhaps up to 
100 GeV. 

The prospect of calculating the baryon asymmetry in 
terms of known experimental parameters make the search 
for a mechanism which amplifies (1) attractive. Recently, 



Farrar and Shaposhnikov Q pointed out one very inter- 
esting mechanism involving very low momentum quasi- 
particles in the quark plasma. In their calculation, the 
scale A is provided by a light quark mass, enormously 
amplifying dcp ■ However it appears that when damping 
due to stronff scattering is included a negligible asymme- 
try results |@,||. 

One might conclude that generating the observed 
asymmetry in the standard model is hopeless. We be- 
lieve this is premature, and in this letter propose a new 
mechanism through which CP violation in the KM ma- 
trix may be greatly amplified in the standard model (or 
extensions such as the minimal supersymmetric version) 
through macroscopic physical effects qualitatively differ- 
ent than those previously considered. 

There is another very large number involved in the 
electroweak transition, namely the horizon Rh in units 
of T~^. The electroweak phase transition proceeds via 
bubble nucleation and growth. Calculations based on 
the perturbative potential show that bubbles fill space 
when their typical radius is ^ 10""* — IO^^Rh- If the 
phase transition is more strongly first order, this could 
increase to a value closer to Rh ~ lO^^T^^. 

Our scenario makes use of this as follows. As bubbles 
form and grow, collisions with fermions in the plasma 
lead to a CP violating instability on the bubble wall, 
namely a condensate of longitudinal Z bosons. The Z 
field can point cither 'out' or 'in' (Figure |l|). Patches 
of each phase initially cover each bubble, separated by 
one dimensional phase boundaries, which are in fact 
Z-magnetic flux tubes. As the bubble expands, small 
patches shrink away under the tension of the phase 
boundary, while larger patches grow with the bubble. 
The long time available before bubbles collide allows each 
bubble to be covered by macroscopic regions of the 'in' 
and the 'out' phases. Even if each bubble becomes com- 
pletely ordered, macroscopic regions of both phases are 
placed in competition when large bubbles collide as the 
transition completes. 

Now as a given phase moves through the plasma, its 
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terminal velocity is given by Vt — where AP is the 
pressure difference between the false and true vacua and 
where the frictional drag/unit area is given by Fjrag = 
Fv. As a result of the CP violation in the KM matrix, we 
expect the 'in' and 'out' phases to feel slightly different 
r's, p'''~p~ ^ dcp- It follows that the two phases move 
through the plasma with slightly different velocities \ v+ — 
V- \ ~ dcp- Now consider the neighbourhood of a phase 
boundary. When the bubbles are large enough, we may 
take them to be planar as far as the local dynamics of 
the phase boundary go. As the bubbles grow, the 'out' 
phase begins to bulge over the 'in' phase, the height of 
the bulge being given by ~ dept. After a time of order 
one expansion time, this becomes comparable to the scale 
which enters in determining the shape of the bubble wall, 
I = where a is the surface tension. The bulge then 
gives rise to a net tangential force on the phase boundary, 
causing it to sweep across the bubble at a velocity of the 
order of the bubble wall velocity. Thus the faster moving 
phase will (again in a time of the order of an expansion 
time) , completely overtake the slower moving phase [ p^ . 
In this scenario, as the electroweak phase transition nears 
completion, bubble surfaces of one CP violating variety 
predominate. Notice that the only place the small CP- 
violating parameter comes in is in ensuring that the 'out' 
and 'in' phases expand at slightly different rates. 
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FIG. 1. Bubbles of the broken symmetry phase growing at 
the electroweak transition, with a condensate of longitudinal 
Z bosons on their surfaces. When large bubbles collide, 'in' 
and 'out' phases are placed in contact, leading to competi- 
tion between macrocopic phases, which greatly amplifies CP 
violating effects. 

The longitudinal Z condensate leads straightforwardly 
to the production of a baryon asymmetry through both 
'local' and 'nonlocal' mechanisms 12|, because 

it couples to fermions in exactly the same way that the 
CP violating phase does in two-Higgs scenarios. 

The remainder of this letter will be devoted to a 
demonstration that longitudinal Z condensate can form 



on propagating bubble walls. The longitudinal Z may 
be defined by the formula i(p*'D^ip = —2gA4>'^ Zfj_, where 
0= VV^Vand gA= j{g1+ 92)^ , with gi,g2 theJ7(l)y 
and SU{2)\Y gauge couphngs. We now consider the clas- 
sical field equation for Z^ , in background of a bubble wall 
propagating through the plasma. In the wall rest frame 
we take the Higgs field to have a fixed profile 4>{z), its 
width Lw being of order the inverse of the Higgs mass 
rriH in the broken phase. The wall is treated as planar, 
which is a good approximation once a bubble gets large. 
We consider a Z field = (0, 0, 0, Z{z)). We may con- 
sistently drop the VF* boson and photon fields, because 
the fermion currents induced by the Z are neutral and so 
only source the Z. Likewise the Z makes no contribution 
to the SU{2) currents in the absence of a W field. 

Next, we drop spatial derivative terms, which are sup- 
pressed by powers of mH/rnz, and iriH/mtop relative to 
the Z mass term and the top quark currents which dom- 
inate. The equation of motion then becomes 



Z = JziZ) + Jf{Z) 



(2) 



where Jz{Z) = —Ag^cjPZ is the current carried by the 
Z condensate and Jf{Z) — gA'4'l'^l^'4' + dVpl^'ip is that 
carried by the fermions. 

In the absence of the bosons, CP invariance holds, 
and so Z — Z is a symmetry Thus there is 

a solution Z ~ Jp = 0. But we shall show that 
this is unstable, so that a nonzero Z current develops. 
We expect this instability to lead to one of two stable 
solutions, related by Z — > —Z, in which Jf{Z) and 
JziZ) are equal and opposite so that the net current 
is zero. The condition for an instability is simply that 
miff = 4.gl(f>^ - idJF/dZ)\z=o < 0. 

The fermion current Jp is calculated by solving the 
time independent Dirac equation in the presence of (jj 
and Z condensates. We have used the WKB approxima- 
tion as a guide to the behaviour expected, and checked 
it against results obtained with an exact linear response 
function. 

The relevant Dirac equation is 



(z7^(9^ + iZ^igv + gAl')) - miz))ij = 



(3) 



where m(z) — y4>{z) with y the Yukawa coupling We 
ignore for the moment the one loop thermal contribu- 
tions to the fermion self energies. From (|) the following 
dispersion relation is found 



= ^pI + (VCpT- gvZr + m2(z) T gAZ{z)r (4) 



for the two eigenstates, which have spin S*^ = ±-i in a 
frame where p± = 0. In the rest frame of the wall, both 
the energy lo and the transverse momentum p± are con- 
served, and the local value of the momentum varies ac- 
cording to (^. It is related to pz^^oo, the value of the mo- 
mentum at z — —00, by the formula E — \/p\ + pI^-od- 
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±i excitations see a 'barrier 




FIG. 2. Profile of a bubble wall, showing the barriers seen 
by particles with spins 5*^ in the positive and negative z direc- 
tions. Positive spin particles incident from the left get further 
than negative spin particles, causing a positive chiral current 
which acts to enhance the Z condensate. This is the origin of 
the instability. 

The fact that left and right handed fermions carry dif- 
ferent charges now comes into play - particles of given 
S'^ carry the same current when travelling in either di- 
rection! One can anticipate the destabilising effect - for 
positive Z there is a region where only > parti- 
cles penetrate from the left, creating a positive chiral 
current. This acts through the equation of motion (|^) 
to further destabilise the Z. Including currents from an- 
tiparticles, which are given by the substitutions J ^ — J; 
gA,v ~^ ~gA,v, multiplies the axial current by two, while 
cancelling the vector current. The parity violation in the 
Z-fermion coupling is crucial - the charged vector bosons 
do not destabilise the Z because the Z/W^ equa- 
tions of motion are parity invariant so that left and right 
moving modes carry opposite currents. 

A more detailed analysis reveals a competing stabilis- 
ing effect which exactly cancels the leading destabilising 
effect at zero wall velocity. The current carried per mode 
is given by the classical formula Jz — ~{duj/dZ) [^ . 
Integrating over momenta this is (for S^- > 0) 



Jf = - 



9 A 



(Pp duj+ 
sfip) 



(2^) 
d^p 



dZ 

dpz 



(27r)3 ^dpz 



-)f{P- 



yz.—oo 



(5) 



(6) 



(the phase space density / is constant along particle 
trajectories). The only Z dependence is now in the Ja- 
cobian, which is just the ratio of the group velocity at 
infinity to that locally. 



duj _ a/ {Pz-oo + gAZy - m? ,Pz -oo 



dpz 



-qaZ 



(^^) (7) 



The Jacobian represents the enhancement of the lo- 
cal particle density due to a 'slowing down' effect. For 
positive Z, S*^ > particles see a 'well' and speed up, 
decreasing the chiral current, and 5^ < particles see 
a 'barrier', slowing down and adding a negative chiral 



current. Thus particles passing 'over' the barrier (in ei- 
ther direction) act to stabilise the Z condensate, whereas 
particles 'bouncing' off the barrier tend to destabihse it. 

Now we compute the net chiral current in the WKB 
approximation. We assume the barrier is monotonic 
{\{dZ / dm)\ < 1), and sum over particles and antipar- 
ticles of both spins, incident from both sides of the wall 
with thermal distributions at z = ±oo. As already men- 
tioned, the leading effect, a square root divergence, can- 
cels at zero - there is a factor of two in the 'bouncing' 
contribution, cancelled by a minus one from particles go- 
ing over the barrier in each direction. But at finite 
more particles are encountered from the left, enhancing 
the number of 'bouncing' particles. The dominant term 
occurs as m(z) approaches moo, and to lowest order in 
Vw is given by 



2gANc ml^Vy, r- — 



+ gAZY-m{zf-{Z^-Z) 



(8) 



where we sum over colors Nn- Note that 



dJp 
dZ 



diverging as z — > oo, so that (within the WKB, free par- 
ticle approximation we have made) the fermionic contri- 
bution to the effective Z squared mass diverges, due to 
the 'bouncing' effect, and an instability always develops 
sufficiently far behind the wall. 

Several comments are in order. The destabilising term 
is proportional to m,^, so the top quark dominates. We 
ignored the decay process t — > W + b, which occurs 
once the t gets a large enough mass on the wall. The 
rate for this process is j^a2ml/mfy{l — m^/mj)^(l -I- 
2m^/mf) ^ nit /no, and the timescale is long compared 
to the time the top quarks spend on the wall, in the 
regime of interest. We have treated the fermions bounc- 
ing off the wall as free particles, ignoring their interac- 
tions with the plasma. More detailed calculations are 
needed to reveal whether these interactions strengthen 
or suppress the destabilising term. One estimate is to as- 
sume the relevant top quarks are emitted with a thermal 
distribution a distance D away from the point z. Then 
in (^) the exponential is replaced by g-^/^^™, with D a 
rather short diffusion length, smaller than (we need 
the exponent to be ~ 3 for instability). But this is too 
pessimistic - it ignores the 'pile- up' of tops in front of the 
wall, particularly pz below the barrier. In any case, a 
full scmiclassical calculation including strong scattering 
in the Boltzmann equation appears quite feasible. 

We emphasise that unlike the calculations of 0| 
and Q which are one dimensional, our calculations are 
three dimensional, and phase space is dominated by large 
perpendicular momenta p± ~ 2T. At large p the quark 
damping rate p3[ is smaller, and in the broken phase 
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the top's large mass makes it less sensitive to collisions 
at low momentum transfer. Note also that a difference 
from Farrar and Shaposhnikov is that the spontaneous 
CP violation here is essentially a classical effect, and 
does not rely on the quantum mechanical coherence. We 
expect it to be less susceptible to destruction by strong 
scattering. 

The WKB calculation is indicative of an interesting 
effect, but since the instability occurs near (and because 
of!) a classical turning point, it is important to check it 
in a full quantum mechanical calculation. We used an 
exact solution to the Dirac equation in the background 
(f>{z)'^ oc (1 + e~°^)~^ in hypergeometric functions, in 
order to compute the full linear response kernel K(x, y) ~ 
{5JF{x)/SZ{y)) Q, which enters the Z equation as 



T^Nc 
47r2 



dyK{z,y)Z{y)) (10) 



The full quantum problem is nonlocal, and the condi- 
tion for instability now reads 



J J dxdyZ{x)K{x,y)Z{y) IQn 



2^2 



JdxZ^x){r{x)/c^l,) 



NrT^ 



100 (11) 



for some trial function Z(x), where we used the bound 
0OO ~ T/g2 ~ 1.5T at the transition - if 4>oo were 
smaller, baryons would not survive in the broken phase. 

We have computed the integrals in ( pT| ) numerically, 
using Mathematica to evaluate the kernel. In order to 
focus on the 'bouncing' effect, we consider the case where 
the only particles incident are those below the barrier 
from the left, and we include both spins, particles and 
antiparticles. We integrate over < < cxd with a 
thermal distribution function and then numerically over 
Pz.-oo, ignoring a phase space factor xLog(l + e~^), with 
X = pz-oo/T which is of order unity, and in any case 
could be enhanced by a chemical potential. We then 
compare the left hand side of ( |ll]) with the leading term 
in the WKB approximation for exactly the same quantity. 
(In WKB, the neglect of particles passing over the barrier 
has the effect of doubling the destabilising term). 

In units where T — 1, we used a — .1 (wall width — 
10), and moo = 1 for the top mass in the broken phase. 
With a Gaussian trial function exp(— (x — Xof /np'), and 
xq = 6Lw, the left hand side of ( [ll|) is greatest for 
w « 2L\Y, and has a value 179, compared to the WKB 
value of 161. The quantum problem thus shows a desta- 
bilising term even greater than the leading WKB approx- 
imation. A comprehensive set of results for different wall 
thicknesses will be presented in p4| . 

Several outstanding questions remain for future work. 
Exactly how docs the CP violation from the KM ma- 
trix enter, and distinguish the 'in' and 'out' phases of the 
Z condensate? In a one dimensional treatment (which 
likely overestimates the effect), one can see from the cal- 
culations of S that a Z independent term ^ dcp enters 



in the equation of motion (g) , and we expect to find sim- 
ilar CP violating terms at other even powers of Z. We 
are investigating the analogous three dimensional calcu- 
lation. Finally, before the baryon asymmetry may be ac- 
curately calculated, we need to find the stable solutions 
for the Z condensate. 
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